Thermodynamics of Strongly Correlated One-Dimensional Bose Gases 
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We investigate the thermodynamics of one-dimensional Bose gases in the strongly correlated 
regime. To this end, we prepare ensembles of independent ID Bose gases in a two-dimensional optical 
lattice and perform high-resolution in situ imaging of the column-integrated density distribution. 
Using an inverse Abel transformation we derive effective one-dimensional line-density profiles and 
compare them to exact theoretical models. The high resolution allows for a direct thermometry of the 
trapped ensembles. The knowledge about the temperature enables us to extract thermodynamic 
equations of state such as the phase-space density, the entropy per particle and the local pair 
correlation function. 
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One-dimensional (ID) systems exhibit a number of pe- 
culiar features which distinguish their physical behaviour 
from those of higher dimensional systems. The precise 
balance between the interaction energy, the kinetic en- 
ergy and the thermal energy strongly affects their proper- 
ties and give rise to a rich landscape of different regimes. 
Theoretically, they can be described for all temperatures 
and interaction strengths by exact models [TJ [5] , solvable 
with powerful numerical methods [3] and fundamental 
mapping theorems [4]. Ultracold atomic gases in one- 
dimensional trapping potentials provide a unique exper- 
imental platform to realize those systems, thus ideally 
complementing this extensive theoretical framework. 

Various physical properties of ID Bose gases have 
been studied experimentally, e.g., pair correlations [5]- 
|S], the momentum distribution [SHU], quench dynamics 
[T2"] and full counting statistics [T3]. However, an ex- 
haustive study of the thermodynamic properties of those 
systems is still missing. This is of fundamental impor- 
tance since the exact knowledge of thermodynamic vari- 
ables and the use of precise thermometry are instrumen- 
tal to study thermalization processes, which in ID are 
known to be strongly reduced |14j . Moreover, thermody- 
namic quantities can be calculated for homogenous sys- 
tems for all interaction strengths [2] allowing for a strin- 
gent parameter-free comparison with the experimental 
results. Until now, these studies have been restricted to 
the weakly interacting regime |15j using magnetic micro- 
traps [T5HT§] . The strongly interacting regime can be 
reached with the help of two-dimensional optical lattices 
[7HH1 HO]- However, in this case only averaged quanti- 
ties have been measured such that the information on 
individual ID gases is smeared out. Only recently, first 
temperature measurements of averaged profiles have been 
reported by our group [8]. Thermodynamic studies of 
strongly interacting ID gases would open new possibili- 
ties to investigate the interplay between non-equilibrium 
dynamics, thermalization, and strong correlations in a 
quantum many-body system. 

In a ID atomic gas (radial trapping frequency u) r ) the 
effective ID interaction strength is given by <7id ~ 2hauj r 



[21] . Here, a is the 3D scattering length. The ratio be- 
tween thermal and kinetic energy |22| is the interaction 
parameter 7, 
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where m is the mass of the particle and niD the 
line-density In the weakly correlated quasi-condensate 
regime the interaction parameter is 7 < 1 while in the 
strongly correlated regime it is 7 > 1. 

Here, we study the thermodynamics of strongly 
correlated ID Bose gases in a two-dimensional (2D) 
optical lattice. Using an inverse Abel transformation 
we extract effective line-density profiles from integrated 
in situ density profiles. We compare the line-density 
profiles to exact solutions and discuss the role of tem- 
perature and interactions. Due to the inhomogeneous 
axial trapping potential u; ax , the temperature is encoded 
in the thermal wings of each profile. This makes an 
in situ thermometry possible and allows us to extract 
thermodynamic equations of state [23 30 . 

We prepare a Bose-Einstein condensate (BEC) of 80 x 
10 3 87 Rb atoms, trapped in a single focused C02-laser 
beam. We use scanning electron microscopy (SEM) 
[311 I32j to manipulate and probe the column-integrated 
density distributions. The SEM emits a focused electron 
beam (EB) (I = 60 11A, E = 6keV, 240(10) nm FWHM) 
which is moved in a rectangular pattern over the atoms. 
To control the number of atoms, we force an evaporation 
by scanning over one wing of the BEC. This addition- 
ally cools the sample to temperatures as low as 10 nK. 
Subsequently, the BEC is loaded within 210 ms into a 
blue-detuned 2D optical lattice with A = 774 nm. 

In order to access different regimes in ID, we pre- 
pared three samples (A,B,C) with initial atom numbers of 
(10(2), 20(3) and 60(5)) xlO 3 atoms. The radial trapping 
frequencies in the optical lattice are u) r = 2tt x 56 kHz (A 
and B) and ui r = 2-7T x 42 kHz (C). Simultaneously to the 
adiabatic ramp of the lattice, we lower the axial trap- 
ping frequency to oj ax = 27r x 8.8 ± 1.3 Hz (A,B) and 
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FIG. 1. (color online), (a)-(c) In situ density profiles from 
data set C (light solid) for increasing distance from the trap 
center, corresponding to different interaction strengths and 
same temperature. The exact YY fits are shown as dark solid 
lines. N denotes the number of atoms. The peak densities are 
n 1D (x = 0): (a,b,c) = (6.2,2.7,0.7)xl0 6 m _1 . The asymptotic 
solutions with the same number of atoms as for the exact YY- 
fit are also shown: TF- profile (T = OnK) (dashed) and ideal 
Bose gas (dotted). 



w ax = 2tt x 12.7 ± 1.7Hz (C). The angle between the lat- 
tice axes and the EB is 45 ° . Therefore, the pixel size is 
chosen to be 273 nm and the total imaging duration is 
30 ms. The fast scanning direction is oriented along the 
ID gases, such that the scan speed is much faster than 
the speed of sound. All single-shot pictures are corrected 
for angle and position drifts and summed up. This yields 
an image sum P for each of the three samples (A,B,C) 
containing 3200, 3900 and 1900 pictures respectively. 

We decompose the integrated density profiles using an 
inverse Abel transformation A~ l . Altough the 2D lattice 
in our setup has a four-fold symmetry it is smeared out 
due to the SEM imaging settings and the post-processing. 
Therefore, the prerequisite of cylindrical symmetry is ap- 
proximately fulfilled. To perform the Abel inversion, we 
make use of the BASEX-method [33] in a modified way 
[34] . The noisy central region [35] (±3 pixel) is inter- 
polated with a Abel inverted gaussian fit on P. Every 
horizontal line in the resulting image R = A~ 1 (P) cor- 
responds to an average of all ID gases which are at the 
same distance from the symmetry axis and thus have the 
same central interaction parameter 70 = j(x — 0). For 
every line in R we perform a fit with the exact Yang- 
Yang theory (YY) [2J, making a local density approxi- 
mation [22J: jjl(x) = /j,o — V ax (a;), where fx(x = 0) = /Xo is 



FIG. 2. (color online), (a) Central interaction strength 70. 
(b) Temperature T from a fit with the exact YY theory. Every 
point represents one line in R for the corresponding data set. 
Negative (positive) values of the x-axis represent the upper 
(lower) half in R. 



the central chemical potential and V ax (x) — muj 2 x x 2 /2. 
The line-density is fixed by a normalization with the 
atom number and the pixel size, leaving the tempera- 
ture as the only free parameter. As can be seen in Figjl] 
the fits reproduce the density profiles very well. The 
temperatures show only moderate variations which are 
compatible with the estimated error (FigJ2|p). We find 
(A,B,C): (T = 11(2) nK, f = 20(4) nK, T = 25(4) nK). 
This indicates an adiabatic loading of the lattice with- 
out significant perturbations. The residual variations of 
T originate from the inversion method as well as the in- 
terpolation in the center. The temperatures were fur- 
ther cross-checked via a fugacity analysis by fitting a 
thermal distribution to the wings of each profile. As 
shown in Figjl] the density profiles change drastically 
with the interaction parameter 79. This is due to the re- 
duction of interaction energy E lnt ~ «id5id- Note, that 
even though 70 is increasing towards the outer tubes, 
the absolute value of the interaction energy drops as 
tt-id- The critical density at which the thermal energy 
dominates is defined via the dimensionless degeneracy 
temperature t(x) — T/T^x) with Td(x) = h 2 n(x) 2 /2m 
[2"2] . For Fig[lji the value in the center is r(0) = 0.2 
and the density profile is close to a Thomas-Fermi dis- 
tribution. For the high-temperature region (r(0) = 15, 
Figjl];) the effect of interaction is masked, because the 
mean inter-particle distance is larger than the thermal 
de-Broglie wavelength and the system is dominated by 
the thermal energy, resulting in a thermal distribution. 
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FIG. 3. (color online). Simplified phase diagram in the 7 — t 
plane for a uniform ID Bose gas with repulsive contact in- 
teraction [55]. The regimes are: (SI) strongly interacting 
regime, (SI high-T) high-temperature regime, dominated by 
thermal energy, (BG) nearly ideal Bose gas and (QC) quasi- 
condensate. The dashed line shows the degeneracy tempera- 
ture r = T/Td = 1. The colored regions show the different 
data sets, whereas the spread in t is given through the tem- 
perature error from the fit. The stars (*) indicate 70 for the 
profiles shown in FigFHand the diamonds (0) for the profiles 
in Fig|4] 
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FIG. 4. (color online). Comparison of selected density profiles 
from all three data sets with same atom number N — 120(10) 
and different temperatures with the theoretical prediction of 
the YY theory at T = OnK. The atoms agglomerate in the 
center for lower temperatures. 



The border (r = 1.2, FigjTb) between the two regimes is 
not sharp and as a consequence, the increasing effect of 
thermal energy towards the wings is smooth. To high- 
light the role of the temperature in the strongly interact- 
ing regime it is convenient to normalize T to a density- 
independent energy scale T\y> — mg 2 D /2/i 2 fcB, where fee 
is Boltzmann's constant [22] • To enter the strongly in- 
teracting regime (SI), the dimensionless temperature t 
has to fulfill t = T/Tid < 1, which is the case for all 
data sets. Using the dimensionless parameters 7 and i, 



we can depict a phase diagram for a uniform ID Bose 
gas [22], as shown in Fig |3] together with the experimen- 
tally accessed regimes. Data set C covers all regimes 
that are below t = 1, from the weakly correlated quasi- 
condensate regime (QC), crossing the SI regime to the 
high-temperature SI regime. To visualize the effect of the 
temperature, we now compare density profiles with dif- 
ferent temperatures but with the same number of atoms 
N = 120(10) (FigJ4j). In order to account for different 
axial trapping frequencies, the axial coordinate x is nor- 

V3, 



For comparison, a YY 



malized by i? Tf = (^) 
profile at T — is shown. With decreasing temperature, 
the density enhances in the center and converges to the 
T = prediction. 

Finally, our method allows for the extraction of ther- 
modynamical quantities. Every pixel in R is an indepen- 
dent measurement of nm(g, n(x),T), yielding access to 
thermodynamic equations of state [23) . The phase space 
density T> — AdB^iD follows as a direct result from R 
(FigJ5^,). Additionally, integrating mo (g, ^{x),T) yields 
the local pressure of the system 



P(giB,^,T) 
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where the integral can be approximated by a discrete 
sum. For a system in local thermal equilibrium, the 
Gibbs-Duhem relation dP = nd[i + sdT is applicable, 
where s is the entropy density. Using eq.pl), the entropy 
per particle S reads: 
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where the derivative is taken for fixed values of <7id and 
/1. To compute the derivative in eq.([3]), we use the lo- 
cal pressures obtained from sets A and B. The derivative 
can be approximated with the finite difference quotient 
(Pa — Pb)/{Ta — Tb) and nxo(n) as the average den- 
sity. The result is shown in FigjSJa together with the 
prediction of the YY theory for an intermediate tem- 
perature of 18 nK. The minimum value of 0.4 x kg is 
comparable to the entropy per particle for a 2D Mott in- 
sulator of 0.3 x ks, reported in [3B|. A lower value of 
0.06 x ks, but for a weakly interacting 2D Bose gas was 
reported in [27 . The agreement with the Yang- Yang the- 
ory is remarkable for large values of ^1, while for smaller 
values, the spread in temperature does not allow for a 
precise determination. Another quantity of interest in a 
strongly interacting system is the local pair-correlation 
function <?( 2 )(0). For many-body systems with repulsive 
contact interaction, the interaction Hamiltonian reads 
Him = a f-J^(x)^(-x)^(x)^(yi)dx, where *(x) and 
, Ff(x) are the field operators. Using the Hellmann- 
Feynman theorem [37] [35] , one can show [55] that g^ (0) 
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FIG. 5. (color online), (a) Phase space density T> for data set A and C together with the exact YY prediction (dashed/solid), 
(b) Entropy per particle S/ks with YY prediction (dashed) and calculated with set A and B (equal <?id)- The shaded area 
shows the error including the temperature uncertainty, (c) Local pair correlation function g^ 2 \0) calculated from sets B and C 
(similar T) with YY prediction (dashed). The shaded area shows the result for different trapping frequencies in the range of 
u ax = 2tt x (12.7±0.05)Hz. 



in a homogeneous system is then given by 



7(%fVi = 
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We approximate eq.Q by using_P(<?iD, fi, T) of set B and 
C. The results are shown in FigJ5p. The measured g^ 2 \0) 
underestimates the predicted result of the YY model by 
50%. This is due to the small relative difference be- 
tween the two pressures of less than 10% which makes 
the procedure very sensitive to little deviations from the 
ideal density profile and amplifies the influence of the ax- 
ial trap frequency and the finite temperature difference. 
Nevertheless, the minimal value of </ 2 )(0) « 0.3 clearly 
indicates strong anti-bunching for a large range of values 
of the chemical potential - as expected for a partially 
fermionized system. 

In conclusion, we have measured effective in situ den- 
sity profiles of strongly correlated ID Bose gases. The 
density profiles show excellent agreement with the Yang- 
Yang thermodynamic theory. The high precision allows 
for a direct determination of the temperature of the 
atoms in the two-dimensional optical lattice. We discuss 
the role of temperature and its interplay with the atomic 
interactions and derive several thermodynamic equations 
of state: the phase space density, the entropy per parti- 
cle and the local pair-correlation function. Our approach 
paves the way to study in situ thermalization processes 
and dynamical properties of ultracold atomic gases in 
optical lattices as well as in bulk systems. 
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